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Several stochastic processes with virtual particles in two dimensional space-time 
are presented whose mean field equations coincide with Schrodinger, Dirac, Klein- 
Gordon and the quantum mechanic equation for a photon. These processes could be 
used to detect discrete space-time features at the Planck scale. 
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I. INTRODUCTION 

With the appearance of interference effects in the early days of qnantum mechan¬ 
ics it became clear that the dynamic of qnantnm systems could not be simulated by 
particles with classical properties such as definite position and momentum, point 
interaction, sharp energy, etc. and the wave representation was preferred. How¬ 
ever with the advent of quantum held theory the particle picture reappeared, not 
as classical entities, but as virtual particles with ephemeral existence, created and 
annihilated at all space-time points and with all energy and momentum values not 
satisfying the usual energy-momentum relations, that is, “ojf the mass shell”. This 
theory allows for an interpretation of quantum mechanics where the virtual particles 
are assigned some ontological reality as components building the quantum £eld[l|. 
In this interpretation, Feynman diagrams are not only a mathematical term of a 
perturbation expansion but represent really existent processes. 

Historically, virtual particles appeared when the solutions of the dynamic equa¬ 
tions of quantum mechanics were expanded in “normal modes” that were interpreted 
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as quantum excitations. In this work we take the opposite way, that is, we assume 
some types of virtual particles with some type of interaction an we show that their 
collective behaviour is described by the well known dynamic equations of quantum 
mechanics. For simplicity, and considering possible computer simulations, these 
“quantum stochastic processes” are studied in two dimensional space-time that in¬ 
clude the essential physical features. Furthermore we concentrate on free systems, 
that is, without external sources. 

The usual quantum mechanic equations are presented as arising from different 
energy-momentum relations, including the less well known quantum mechanic equa¬ 
tion for the photon. These are specialized for two dimensional space-time and several 
stochastic processes are presented whose mean held equations coincide with quan¬ 
tum mechanic ones. Finally some speculations concerning possible interesting non 
equilibrium effects of the processes are suggested. 

II. QUANTUM DYNAMICS FOR FREE SYSTEMS 

The energy-momentum relation of a free physical system determines its dynamics 
because energy and momentum are the generators of space and time translations. 
For nonrelativistic systems the energy-momentum relation is 

2m 

where E is the energy and = P ■ P is the momentum amplitude squared. In the 
relativistic case the relation is (in units such that c = 1 and with signature -I-) 

. ( 2 ) 

Some inconveniences due to the quadratic nature of the expression above are avoided 
by postulating a linear energy-momentum relation, that must also satisfy the rela¬ 
tivistic quadratic relation. This can only be achieved with coefficients 7 ^ satisfying 
a Clifford algebra { 7 ^, 7 ^} = 2r]'^''. Then we have 


— 7 ■ P = m . 


( 3 ) 
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The relativistic relation (EJ is also valid for massless particles, m = 0, and in this 
case the energy becomes a linear dependence with the momentum amplitude 

^ = |P| . (4) 

For massive particles, spin S and momentum P are decoupled and can have arbitrary 
directions; however, for particles with zero mass, like the photon, special relativity 
requires that spin and momentum must be collinear (transversality constraint), that 
is, 

S X P = 0 , (5) 

and the sign of S ■ P corresponds to the positive or negative helicity of the photon. 
The absolute value in Eq. (jl]) is problematic when we transform the observables into 
operators. In order to avoid the absolute value, it is convenient to combine the 
linear dependence of energy and momentum with the two possible helicity states of 
the photon by choosing the energy-momentum relation for the photon as[^ 

E = S-P . (6) 

When this quantity is negative, it must be interpreted as a photon with positive 
energy but with negative helicity. 

In quantum physics, the energy and momentum observables are represented by 
operators in a Hilbert space and when we choose for it the space of square integrable 
functions, where E —)■ idt and P —)■ —iV (in units such that h = 1), we obtain the 
well known Schrodinger, Klein-Gordon and Dirac equations. In the representation 
to. =p.n 1 whe.e the epi.t ope.ato.e S, have mat.ix elements (S,),. = -^e„, the 
quantum mechanics equation for the photon, obtained from Eq.(j6]), is[2[ 

idtEj = SjikdiEk , (7) 

where are the tree components of the photon state in the tree dimensional Hilbert 
space for spin 1. (However, the photon states are restricted to a two dimensional 
subspace spanned by the helicity ±1 states.) 

Notice that one can formally derive Maxwell’s equations in vacuum (that is, 
with vanishing charge and currents) from Eq.([7]); however their interpretation is 
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quite different: the equation above describes the dynamics of a one photon system 
whereas Maxwell’s equations are the evolution equations for a set of observables 
-the electromagnetic field- of a system made by an indefinite number of photons. 
Furthermore the alleged derivation is based on the conceptual error of identifying the 
tree dimensional Hilbert space of spin with the tree dimensional physical space j^. 

III. QUANTUM DYNAMICS IN TWO DIMENSIONAL SPACE-TIME 

For Schrodinger and Klein-Gordon equations, the specialization for two dimen¬ 
sional space-time {t, x) is straightforward: we just make the association E idt 
and P —idx and we get 

dti^ix^t) = , ( 8 ) 

(dt - dl)'4){x, t) = -m^'ilj{x, t) . (9) 

For Dirac equation we choose the two dimensional matrices and spinor 



and we become from Eq. fl^ 

dt'<p2{x,t) = -dx'il)2{xp) - inii)i{x,t) 

dt'ilJi{x,t) = dx'ipi{x,t) - im'ii)2{,x,t) . ( 11 ) 

The quantum equation for the photon in two dimensional space-time is not trivial 
because the transversality constraint ([5]) and equations ([6]) and ([7]) are essentially in 
tree dimensional space. We can however make some formal modihcations of Eq. ([7|) 
in order to get a two dimensional equation. Let us fix the index I at the value I = 1 
and associate di = dx- The indices j, k cant take the values 2, 3 and 3, 2 resulting 
in two coupled equations: 


idt'1'2 = - 4*^3 , 


idt^3 = +dx^2 ■ 


( 12 ) 
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Although the identification of Eq.(I7]) with Maxwell’s equations is questionable, 
but formally possible, in the two dimensional case we can also obtain the two di¬ 
mensional Maxwell’s equations in vacuum with the misidentification of the photon 
states with the electric and magnetic fields: ^2 ^ E and —?■ B obtaining 

dtE = d,B , 

dtB = d^E . (13) 

These equations are also obtained by generalizing the four dimensional Maxwell’s 
equations to arbitrary dimension with the techniques of differential geometry and 
then specializing it to two dimensional space-time|3|. 

The set of equations flT^ can be decoupled with the replacement Er = E 2 — 
and El = E 2 + obtaining 


QiEr — —OxEr , 

dtEL = +dxEL . (14) 

If we try solutions with the form f{x — vt) we can see that Er and Er represent 
a photon moving to the right {v = 1) or to the left {v = —1) and therefore they 
are the one dimensional counterparts of the tree dimensional right or left handed 
photons. Notice that the H— signs in Eq. flT^ correspond to different helicity states 
and in Eq. fimi they correspond to different directions of propagation. 


IV. SCHRODINGER PROCESS 

Let us assume two types of virtual particles A and i? in a one dimensional lattice 
with sites at distance A . These particles have an asymmetric interaction (violating 
the action-reaction principle): A particles reject B from its location whereas B 
particles attract A to its location. Metaphorically we can say that A hates B and B 
loves A or that B particles are good and A are evil. The rejection of a particle A or 
B from a site is formalized by the creation of an antiparticle A or S at the site. Each 
site can be occupied by any number of particles of type A i? or by their corresponding 
antiparticles A B. Particles and antiparticles of the same type annihilate in each 
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site of the lattice leaving only the remaining excess of particles or antiparticles of 
both types A and B. At each discretized time step, t —)■ t + 1, corresponding to a 
time evolution by a small amount r, each particle of type A creates two antiparticles 
B in the same site and one particle B in each of the two neighbouring sites. This 
is equivalent to moving one particle B to the right and another to the left. In a 
similar way, B particles move neighbouring A particles attracting them to its site. 
This process is realized with some amplitude (probability, if normalized) k. 

A (B) (2B) (B) 

B (A) (2A) (A) . (15) 

The same reactions occur exchanging particles and antiparticles. 

Before we write the master equation for the time evolution, we can notice some 
qualitative features of the process. If A rejects B but B attracts A then, transitively, 
A reject themselves and diffuse. Same conclusion is reached for B. It is obvious 
that the process has diffusion because in each time step particles and antiparticles 
are created in neighbouring sites. For instance, if we start with two particles, one A 
and one B, in one site, after three time steps we will have 110 particles occupying 
seven sites. It is less obvious that, even though the process has left-right symmetry, 
we may also have drift to the right or to the left. In order to see how this is 
possible, recall that A particles reject B particles and B attraet neighbouring A 
particles. Therefore, if we have an asymmetric conhguration like AB, the center of 
the combined distribution will move towards B. The drift direction and velocity is 
then encoded in the shape and relative distribution of both types of particles. It 
is remarkable that, although the distribution of particles are widely distorted after 
few time steps, the drift direction and velocity remain invariant (in the continuous 
limit). 

Let as(t) and bgit) be the number of particles of type A and B respectively at the 
site s at time t. When as{t) or bs(t) take negative values they denote the number 
of antiparticles. At a particular site of the lattice, the number of particles change 
as particles or antiparticles are created there by the particles in neighbouring sites. 
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The time evolution of the process is then dehned by the master equations 
as{t + 1) = as{t) - k (bs-iit) - 2bs{t) + bs+i{t)) 

bs{t + 1) = bs{t) + k {as-i{t) - 2as{t) + as+i{t)) . (16) 


Choosing time and space scales such that ^ = 1 we can write these eqnations as 


asjt + 1 ) - asjt) _ bs-i{t) - 2bs{t) + bs+i{t) 


T 


A2 


bs{t+ 1) - bs{t) _ as-ijt) - 2as(t) + a^+i(t) 


(17) 


T A2 

In these eqnations one can easily recognize the discrete version of the time derivative 
and of the second space derivative. Taking then the limit r —)■ 0 and A —?• 0 with 
^ 1 and k = ^ and also replacing as{t) and bs{t) by continnons fnnctions a{x, t) 

and b{x,t), the eqnations above become 


dta(x,t) = - dlb(x,t) 

2m 

dtb{x,t) = ^ dla{x,t) . (18) 

Now we can combine the real amplitndes of the virtnal particles in one complex held 
'ip{x,t) = a{x,t) + ib{x,t) and both eqnations above resnlt in Schrodinger eqnation 
given in Eq.([H]). 

The virtnal particles process presented here is a mnch simplihed version of a 
similar process stndied in more detaillsjd]. 


V. DIRAC PROCESS 

Let ns assnme now fonr types of virtnal particles A, B,C, D in a. one dimensional 
space. These particles always move and create and annihilate each other. Precisely, 
A and B always move to the left whereas C and D move to the right; D particles 
create A bnt A annihilates D and similarly, B particles create C bnt C annihilates 
B. 

Notice that with regard of the movement, the pair {AB) behaves as chirality 
opposed to the pair {CD), that is, we may think of {AB) as having negative helicity 
(left handed) and {CD) with positive helicity (right handed). Fnrthermore, with 
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respect with creation and annihilation, the pair {AC) behaves as the antiparticles of 
the pair {BD). Therefore the process globally contains the essence of antimatter but 
not in an identihed way, that is, none of the particles is individually the antiparticle 
of another. Another interesting qualitative feature is that since D creates A but A 
annihilates D, transitively, D has self annihilation and the same can be said for all 
other particles. This self annihilation component is essential to the virtuality of the 
particle: they can not have permanent existence. 

In order to formalize this process let us assume again an inhnite one dimensional 
lattice with sites at distance A occupied by some number of these particles. Each site 
can be occupied by any number of particles of type A, B, C, D. At each discretized 
time step, f —)■ f+1, corresponding to a time evolution by a small amount r, particles 
move one lattice site and are created and destroyed with some amplitude k. In this 
case we choose space-time units and scale such that y = 1- 

Let as{t), bs{t), Cs{t), ds{t) be the number of particles of type A, B, C, D 
respectively at the site s at time t. At a particular site of the lattice, the number of 
particles change as particles migrate, are created or destroyed. The time evolution 
of the process is then defined by the equations 


as{t + l) = as+i{t) + kds{t) 

bs{t + 1) = bs+i{t) - kcs{t) 

Cs(t-hl) = Cs-i{t) + khs{t) 

ds{t -h 1) = ds-i{t) - kas{t) 


(19) 


We will see that this simple process is described by Dirac equation. Let us first 
subtract the same quantity from each side of these equations and then use the 
relation y = 1 to obtain 


a^jt -H) - as{t) 

T 

bs{t -H) - bs{t) 
T 

Cs{t -H) - Cs{t) 

T 

ds{t -hi) - ds{t) 


®s+l(^) ^s{^) 

A 

bs+i{t) - bs{t) 

A 


+ -^ds{t) 

k , , 

- jCs{t) 


Cs{t) - Cs-l{t) 

A 

ds{t) - ds-i{t) 

A 


+ -j^bs{t) 

k , , 

- jas{t) . 


T 


( 20 ) 
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In these master equations one can easily recognize the discrete version of the 
space and time derivatives. Taking then the limit r—)-0, A^O, A:—)-0 
with ^ ^ 1 , I —)■ m and replacing as{t),bs{t),Cs{t),ds{t) by continuous functions 
a(x, t), b{x, t), c(x, t), d{x, t), the equations above become 

dta{x,t) = dxa{x,t) + md{x,t) 
dtb{x,t) = dxb{x, t) — m c{x, t) 
dtc{x,t) = —dxc{x,t) + mb{x,t) 

dtd{x,t) = —dxd{x,t) — ma{x,t) . (21) 

We can now combine the first two real equations in a complex one by adding the 
first with the second multiplied by i (and similarly for the last two) 

dt {a{x, t) + i b{x, t)) = {a{x, t) + i b{x, t)) — im (c(x, t) + i d{x, t)) 

dt{c{x,t) + id{x,t)) = —dx{c{x,t)+id{x,t)))—im{a{x,t)+ib{x,t)) ,(22) 

and in terms of the complex fields t) = t) + ib{x, t) and ' 02 (a^, t) = c(x, t) + 
id{x,t), the expression above results in Dirac equation given in fITT]) . 

Notice that "01 and 02 have opposed helicity and represent the antimatter of each 
other. 


VI. KLEIN-GORDON PROCESS 

After giving two stochastic processes with mean field equations corresponding 
with Schrodinger and Dirac equations we present now a virtual particle process that 
results in Klein-Gordon equation ([9]) arising from the relativistic energy-momentum 
relation (0. 

The second space derivative is obtained in a similar way as in the Schrodinger 
process in Eq. flT^ (but with just one type of diffusing particle) and the rest energy 
term is also simple, but the second time derivative presents some difficulty because, 
in the discrete time version, it involves tree times: f -|- 1, t, t — 1. Notice that 
the Schrodinger and Dirac processes are markovian because the state at time f -|- 1, 
given by the master equations Eqs. fll6lll9p . are uniquely determined by the state at 
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time t, without memory on past times. In the present case we must abandon the 
markovian property because the master equation will also involve the state at time 
t-1. 

Let as{t) denote the number of particles (antiparticles, if negative) of type A at 
site s at time t in a one dimensional lattice with lattice constant A. At each step r of 
discretized time t —)■ t + 1 all particles reproduce and migrate to neighbouring sites: 
oAo Ao A and are annihilated with probability k. The non-markovian property 
is that all particles present at time t — 1 are removed. The master equation is then 

as{t + 1) = as+i{t) + as-i{t) - kas{t) - as{t - 1) . (23) 

2 

With some trivial manipulations, and using ^ = 1, we can write this equation as 

as(f - 1) - 2 as(f)a5 (t-Fl) as-i{t) - 2as{t) + as+i{t) k 
r2 “ A2 A2® 

Now, taking the limit r^-0. A—?-0, fc—?-0 with ^ 1 , ^ rm? and replacing 

as{t) by the continuous held '0(x,t) we obtain Klein-Gordon equation ([9]). 

It is possible to recover the markovian property for a Klein-Gordon process if we 
make an additional postulate. A process described by agif) can be said to be causal 
or markovian if the state at t -|- 1 is uniquely determined by the state at time f, that 
is, it satishes an evolution equation of the type as{t -|- 1) = F[ar{t)]. Now we can 
dehne the process to have reversed causality or reverse markovian if the state at time 
t — 1, that is, not at future but at past time, is uniquely determined by the state at 
time t. Such a process satishes an equation of the type Ogit — 1) = F[ar{t)]. Notice 
that the requirement of reversed causality is different from time reversal symmetry 
T : t ^ —t. Gausality means that the future is uniquely determined by the present 
whereas in reverse causality the past is uniquely determined by the present. Let us 
consider then a markovian Klein-Gordon process described by 

ag{t -h 1) = ^ag+i{t) + ^as-i(t) - k^ag{t) . (25) 

This is a process similar to the one of Eq. but without the non-markovian term 
ag{t — 1). Now, if we require reversed causality we also have 

ag{t - 1) = l:ag+i{t) + ^ag_i{t) - kl-ag{t) . 


(t) (24) 


(26) 
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Adding the last two equations we get 

as{t + 1) + as{t - 1) = as+i{t) + as-i{t) - kas{t) . (27) 

This is the same master equation fl2^ for the non-markovian process that in the 
continuous limit results in Klein-Gordon equation. 

Here we have taken a real held but the same could be done with a complex or a 
many component held in order to account for particles with charges. 


VII. PHOTON PROCESS 

The photon process is the simplest of all. Let us assume two types of particles 
R and L. In each time step t —?• f + 1 of discretized time, corresponding to a time 
evolution by a small amount r, R particles move one site to the right and L particles 
move to the left in a one dimensional lattice with lattice constant A. The particles 
can not be at rest and move with constant speed Let rs{t) and ls{t) denote the 
number of right and left moving particles at site s at time t. Then, the master 
equations are simply 


r,{t + l) = r^_i(t) 

hit + 1) = h+lit) ■ (28) 


Subtracting Vsit) [/s(f)] on both sides of the hrst [second] equation and choosing 
time and space scales such that ^ = 1, we can write these equations as 


rsjt + 1) - rsjt) _ rsjt) - r^-i(t) 
T A 

tsit + ^) — tsit) _ h+lit) — hit) 

T A 


(29) 


Here we hnd the discrete version of time and space derivatives. Taking then the 
limit r —)■ 0 and A —)■ 0 with ^ ^ 1 and replacing rsit) and hit) by continuous 
functions 'R+iixR) and 'RiixR) we obtain the quantum mechanic dynamic equation 
for a photon given in Eq. flTTD . 
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VIII. CONCLUSION 

In a way consistent with the quantum held theory interpretation of quantum 
mechanics we have found stochastic processes with virtual particles whose mean 
held equations correspond to Schrodinger, Dirac, and Klein-Gordon equations, as 
well as the corresponding quantum mechanic equation for the photon. 

One possible interest of studying these “quantum stochastic processes” appears if 
we consider that the mean held equations are obtained when we take the limit from 
a discrete to a continuous space-time. However it is well known in non-equilibrium 
statistical mechanics that there are some features of discrete space-time processes 
that could appear in a computer simulation that are not present in the solutions of 
the mean held equations. This could be interesting in order to predict some ehects 
due to an inherent discreetness of space-time of the order of Planck time Tp 
and length Lp k. as is suggested in some quantum gravity theories. Such 

simulations with t k. Lp and A ~ Lp are beyond reach with today computing power, 
however with the exponentially increasing computing capabilities it may sometime 
be possible. 
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